Observing Collapse in Colliding Two Dipolar Bose-Einstein Condensates 
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We study the collision of two Bose-Einstein condensates with pure dipolar interaction. A sta- 
tionary pure dipolar condensate is known to be stable when the atom number is below a critical 
value. However, collapse can occur during the collision between two condensates due to local density 
fluctuations even if the total atom number is only a fraction of the critical value. Using full three- 
dimensional numerical simulations, we observe the collapse induced by local density fluctuations. 
For the purpose of future experiments, we present the time dependence of the density distribution, 
energy per particle, and the maximal density of the condensate. We also discuss the collapse time 
as function of the relative phase between the two condensates. 

PACS numbers: 03.75.Hh,03.75.Kk 



I. INTRODUCTION 

Since the experimental observation of 52 Cr Bose- 
Einstein condensate (BEC) [l|, there has been a growing 
interest in the study of ultracold dipolar gases. 52 Cr has 
a magnetic dipole moment of /i = 6/is (/ib is the Bohr 
magneton) which has at least 36 times larger dipolar in- 
teraction strength than its alkaline counterparts. There- 
fore, 52 Cr is an ideal choice for investigating novel dipolar 
effects in neutral atoms. It has been shown theoretically 
and experimentally that there are detectable modifica- 
tions to the condensate dens ity pro file [j, 3, HI, IE H, B, B, HI 
and elementary excitations [Iol.llll [T^.ll3l[l4| due to this 
long ranged and anisotropic interaction. 

The stability of a dipolar condensate is a fundamental 
question to answer since the dipolar interaction is par- 
tially attractive. One feature of the dipolar condensate 
is that the effective dipolar interaction depends on the 
shape of the trap. This can be roughly understood from 
a simple argument based on the competition between the 
potential energy per particle u p and the dipolar interac- 
tion energy per particle u c i- Suppose the dipolar gas is 
polarized along the z-axis and confined in a cylindrically 
symmetric trap with an aspect ratio A = u) z /u>x- With- 
out the trapping potential, the dipoles tend to arrange 
into a head-to-tail configuration and lowers Ud which re- 
sults in an unstable condensate. This is also true for a 
prolate trap (A <C 1) because this configuration also low- 
ers Up. However, in an oblate trap (A 3> 1), there is a 
competition between u p and Ud- Although u p is almost 
independent of the atom number N, the magnitude of 
Ud increases as N increases in general. When u p domi- 
nates Ud, i.e. N is below a critical value N c , the dipoles 
are prone to arrange into a head-to-head configuration 
which is stable. In the opposite case, the gas still favors 
a head-to-tail configuration which is again unstable. 

The dependence of this stability on the trap aspect 
ratio and atom number is investigated more thoroughly 
in a recent publication where the authors show that the 
stability diagram exhibits richer physics beyond our in- 
tuitive understanding Q ■ As shown in the stability dia- 
gram from Ref. Q , an increase in A will tend to stabilize 



a dipolar condensate and for a given A, there is always a 
critical value of N c above which the condensate is unsta- 
ble, in agreement with our simple analysis. However, the 
dipolar interaction can cause the formation of a struc- 
tured condensate, e.g. a biconcave one, in addition to 
a normal Gaussian shaped condensate. The underlying 
mechanisms of the instability can be analyzed from the 
Bogoliubov-de Gennes equation which are referred to as 
angular- and radial-roton instability, respectively @, fll| . 

Previous studies concerning collapse in dipolar gases 
have focused on the response of a stationary conden- 
sate to a modulation of the s-wave scattering length 
[1, [H, [H, [It]]- In this paper, we want to investigate 
the possibility of observing collapse induced by purely 
dipolar interaction in a dynamic collision process. For 
this purpose, we will study the collision dynamics of two 
dipolar condensates and discuss the collapse effect. A 
similar scenario of overlapping several independent con- 
densates is briefly discussed where the effect of collapse 
on the interference fringes is observed Q. Another rel- 
evant scenario is for the case of pure attractive s-wave 
scattering where colliding two bright solitary waves may 
also lead to collapse [18| . We will give a more thorough 
study and present more detailed results such as density 
distributions and the collapse time to support future ex- 
periments. The structure of this paper is as follows. In 
Sec. II, we start with the generalized Gross-Pitaevskii 
equation, giving our simulation parameters and the nu- 
merical scheme. In Sec. Ill, we present numerical results 
and talk about the collapse effect in detail. Finally, we 
give the conclusion in Sec. IV. 



II. THEORY 

The dynamics of the two BECs at sufficiently low 
temperature are described by the generalized Gross- 
Pitaevskii equation (GPE), 
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with various terms listed below 



III. NUMERICAL RESULTS 
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where the wave function is always normalized to unity. 
m is the atom mass, N is the atom number, and clq is 
the s-wave scattering length. The trap potential assumes 
the following form 

Vtm P (r,t) = ^mujl(p 2 + \ 2 z 2 )+Ae-^0(-t), 

where p — (x, y) and #(•) is the Heaviside step function. 
It describes a trap potential which at t = is a combina- 
tion of a cylindrically harmonic trap (A = lu z /lu±) plus 
a central Gaussian barrier (with height A and width w) 
along the x-axis and the barrier is removed immediately 
after t — 0. The dipolar interaction potential for a gas 
polarized along the z-axis is 



Vdd(r) = c d - 



3z 2 



Cd = / i oA* 2 /(47r) where fio = 4ir x 10~ 7 T-m/A is the vac- 
uum permeability and /i is the magnetic dipole moment 

of the atom. 

We adopt the length scale ciho = \A/ (muj±) and the 
time scale l/u>± as those of a harmonic oscillator and 
substitute r — > r/ai lo , t — » ta;j_, .4 — > A/(hu±), it) — > 
w/a,ho, Q = Na a /a ho , and Z? = Nc d /{hu>^a\ ) into Eq. 
([T]) . We then arrive at the dimensionless version of the 
generalized GPE 



. dj>(r,t) 
+D [ dr' 



-iv 2 + T/ trap (r,t) + 47rQ|VXr,i)| 2 



n'|5 



3(z - z') 



A2 



r'|5 



|^(r',t)| 2 U(r,i), (2) 



where Vk ap (r,i) = ±(p 2 + A 2 z 2 ) + Ae ^?6(-t). In this 
paper, we study the dynamics of BECs with purely dipo- 
lar interaction, assuming the s-wave scattering length can 
be tuned to zero (Q = 0) by modulating the magnetic 
field. The interaction term (the last term in Eq. ©) 
can be conveniently computed by making use of the Fast 
Fourier transform T with the help of the following iden- 
tities 



F[H d ]= NT[V dd {r)]F[\^{r,t)\ 2 ], 
F[V dd (r)] =c d — (3cos 2 k -l), 



(3) 



where 9^ is the angle between the conjugate momentum 
k and the z-axis. The H d term can be obtained via an 
inverse Fast Fourier transform. 



In this section, we present our numerical results for 
colliding two BECs. We first briefly discuss the initial 
state. Then we will investigate the collision dynamics in 
detail. 

We use 52 Cr atom with /i — 6/is in our numerical 
simulations. For the harmonic trap, the transverse fre- 
quency is fixed at lo±_ = (27r)25 Hz. The correspond- 
ing length scale is ciho = 2.78/im and the time scale is 
= 6.37ms. The mesh in all three directions is 
192 x 0.1. The initial state at t = is obtained first 
by imaginary time relaxation. We then add a phase fac- 
tor e 1 ^ to the right side of the wave function to account 
for a possibility of uncertainty in the initial phases be- 
tween the two BECs [lj| . To verify the reliability of our 
numerical results, we search for the critical atom num- 
bers for the ground state in harmonic traps with different 
aspect ratios and compare them with those reported in 
Ref. @. We find our results for the ground state are in 
good agreement with Ref. Q. For example, we found for 
A = 2 (ut z — (27r)50Hz), the critical atom number for the 
ground state is around iV c = 4400. We also search the 
critical atom numbers for other excited states. For exam- 
ple, for p-wave soliton, the critical atom number is found 
to be larger than the ground state for the same aspect 
ratio. Here the p-wave soliton refers to the state with 
a node along one direction (say x) so the ansatz func- 
tion for the imaginary time propagation takes the form 
of xe~ p / 2 ~ Xz I 2 . The detailed study of excited states 
will be reported elsewhere. In this paper, we are inter- 
ested in the situation of colliding two ground state con- 
densates with a strong confinement along the z-axis so 
we choose A = 10 (u z = (2tt)250Hz) and N = 19000 
where the ground state is found to be stable. We will use 
them throughout our simulation for the collision dynam- 
ics. The column density p(x, y) along y = is shown as 
the red solid curve in Fig. [T]where p{x,y) = J dz\ip(r)\ 2 . 
When a central barrier is added, it is shown as the blue 
dashed curve in Fig. [1] The dipolar interaction strength 
is computed to be D = 16.5. The parameters for the 
central barrier are ^4 = 15 and w = 1 which are chosen 
so that, at the initial time t = 0, there is no substantial 
overlap between the two condensates. We want to em- 
phasize that this atom number is far below the critical 
value N c ~ 70000 estimated from Q for the same single 
harmonic trap, as the motivation is to observe the local 
collapse induced by density fluctuations rather than the 
global collapse induced by an overall critical atom num- 
ber. After t > 0, the condensates are out of equilibrium 
and start to collide with each other and oscillate in the 
trap. We explore the subsequent dynamics using real 
time propagation. 

We first discuss the simplest case of two non- 
interacting BECs (Q = and D = 0), assuming the 
dipolar interaction can be tuned to zero [20j. In Fig. [21 
we show our numerical results of linear density p(x, t) as 
function of time t for the relative phase = 0. p(x, t) 
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FIG. 1: Column density p(x,y) along y = as function of 
x. The red solid (blue dashed) curve is for ground state in 
the single harmonic trap without (with) a central gaussian 
barrier. 



is obtained by integrating along both the y- and z-axis, 
i.e. p(x,t) = J J dydz\ip(r,t)\ 2 . At the initial time, the 
two BECs are well separated. After the barrier is turned 
off, the two BECs start to collide with each other. As 
a consequence, we can see an interference pattern be- 
tween them. The interference pattern disappears when 
the two BECs pass through each other. The two BECs 
are then reflected by the harmonic trap and ready for 
the collision in the next cycle. Such a necklace pattern 
of p(x, t) is expected to keep repeating itself, with only 
a few cycles are selectively shown in Fig. [3] Because the 
two BECs are non-interacting, this pattern persists and 
collapse will not take place. 




FIG. 2: Linear density distribution p(x) as function of time 
t for two non-interacting condensates with = 0. Blue (red) 
stands for the minimum (maximum) density. 

Now we present the numerical results for dipolar con- 
densates (Q — and D = 16.5). Firstly, the density 



distribution differs significantly from that in the non- 
interacting case. Secondly and more importantly, we ob- 
serve the collapse of the two BECs which is absent in 
the non-interacting case. Our numerical results for rel- 
ative phase <fi — are presented in Figs. [3][5] In Fig. 
[3l we show the density distributions p(x,t) and p(x,y,t) 
as function of time t. From the upper figure in Fig. [3l 
we can see that the two dipolar BECs behave similarly 
to non-interacting BECs: they collide with each other, 
interfere, and are reflected by the harmonic trap. How- 
ever, the differences between them still merit some dis- 
cussions here. One difference is that the period of the 
first cycle is about 4.5 which is larger than that of non- 
interacting BECs (~ 3). This slowed motion is reminis- 
cent of the damping effect for non-dipolar condensates 
(Q ^ and D = 0) where the interaction causes com- 
plex flow patterns acting as a damping force [2ll . I22I |23| . 
Here we also find complex flow patterns for the dipo- 
lar interaction. Selected plots for the probability current 
,/ = ^lm[ip*(x,y,z,t)V^(x,y,z,t)] ("Im" denotes the 
imaginary part) are shown in Fig. [3] The upper row of 
Fig. 2] is for non-interacting condensates where the flow 
lines are all along the x-axis. However, in the lower row 
of Fig. H] as for interacting condensates, we can see flow 
lines bent by the interaction. As a result of the damped 
motion, roughly 5 interference fringes can be seen, cf only 
3 interference fringes can be seen for the non-interacting 
case. Another difference is the increased density in the 
central region at the final time (t ~ 6.5) which leads 
to the collapse effect. Since the total atom number is 
much lower than the critical atom number, this collapse 
is purely seeded by local density fluctuations. In this 
case, the interference is responsible for the enhancement 
in the density distribution. Although the dipolar inter- 
action is cylindrically symmetric, the existence of the ini- 
tial central barrier breaks this symmetry. As a result, the 
density distribution will develop an anisotropic pattern 
in the transverse plane. This cannot be seen from the 
linear density but instead can be seen from the column 
density p(x,y,t). In the lower part of Fig. [31 we show 
6 time snapshots (a-f) for p(x,y,t) whose times are also 
marked on the horizontal axis of the upper figure. The 
density patterns in Figs. 3(a)-(d) clearly show a cycle of 
collision. We can see that, due to the dipolar interaction, 
the density pattern of Fig. 3(d) is not the same as that of 
Fig. 3(a). At a later time t ~ 5.6 (Fig. 3(e)), the density 
also exhibits a distribution along the y-axis: the density 
is maximal in the central region and surrounded by 4 
lobes. At the final time t ~ 6.5 (Fig. 3(f)) just before 
the collapse, the central density maximum evolves into a 
singularity and triggers the collapse. Note that this sin- 
gularity is purely artificial since CPE cannot handle the 
post-collapse dynamics. However, the CPE can still pro- 
vide an accurate prediction of the onset of this collapse. 
Therefore, we still choose to present Fig. 3(f) to demon- 
strate the singular density profile. From the definition of 
D, where N/a\ a is just the averaged density, it is perhaps 
not surprising that the local density fluctuations may in- 
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duce collapse even though the total atom number is far 
below the critical atom number. What is somewhat more 
interesting is that the collapse does not happen during 
the first cycle of collision (t ~ — 4.5) but only at a later 
time (t ~ 6.5). This is different from the previous study 
for the case of attractive s-wave scattering where the col- 
lapse happens during the first cycle of collision [l^]. In 
our case, it takes longer time for the anisotropic dipolar 
interaction to build up the local density fluctuation and 
eventually singularity. 




FIG. 3: Upper figure: linear density distribution p(x) as func- 
tion of time t for <f> = 0. Lower figure: column density p(x, y, t) 
as 6 different times labelled by (a)-(f). Field of view in each 
subplot is (x,y) — [—3.5,3.5] x [—3.5,3.5]. In both figures, 
blue (red) stands for the minimum (maximum) density. How- 
ever, the colormap of the upper figure is different from that 
of the lower figure and is adjusted for better visual effect. 

To track this collapse more quantitively, we show the 
energy per particle E and the maximal density p max as 
function of time t in Fig. /? max = max{ r } IVK 1- )! 2 
at each time step. The heating effect arising from 
the sudden removal of the barrier potential is negligi- 
ble as the energy variation at the initial stage is about 
0.5Tiuj±/kB ~ nK, so the description using GPE is still a 
good approximation. It is clear from Fig. [5] that the col- 
lapse is signaled by a sudden drop (increase) in E (p max ) 
due to the singularity developed in the condensate den- 
sity. Here we are only interested in the onset of this 
collapse. Although the collapse dynamics afterwards are 
also interesting, a similar scenario of collapse has already 
been experimentally observed and discussed in detail in 
Ref. [1] and will not be discussed in this paper. 

For other choices of </>, we also observe the collapse 
effect. For example, in Fig. [51 we show p(x, t) as function 
of time t for (j> — it. We find that the period of the first 
cycle is close to that of <j> = 0. Since the condensates 



FIG. 4: Probability current J for (j> = at t = 1.43 (left) 
and t = 2.86 (right). The upper (low) row is for the non- 
interacting (dipolar) case. Field of view in each subplot is 



{x,y) 



-3,3] x [-3,3]. 
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FIG. 5: Energy per particle E (red solid) and the maximal 
density p ma x (blue dashed) as function of time t <j> — 0. 



are always in a phase of it, the two BECs never pass 
through each other. It seems as if they just collide and are 
bounced back from each other. However, we still observe 
strong density fluctuations and the induced collapse at 
t ~ 8.7. In this case, the density maximum is not in the 
center, rather it is found at two different locations at the 
same time. The complete dependence of the collapse time 
t c on <p is shown in Fig. [7] t c is defined as the time when 
|1 - E(t c )t c / fa E{t)dt\ = 8 and we choose S = 5% here. 
Although a purely empirical definition, t c does capture 
the onset of the energy variation due to the diverging 
density profile. We find that a moderate change in 5 
does not change our conclusion qualitatively. The overall 
trend of the curve is that t c increases as <f> increases. This 
can be explained as follows. For small <f> ~ 0, the collision 
behaves as "attractive" and the maximal density is likely 
to be found in the trap center resulting in a enhanced 



5 



density. While for large <f> ~ w, the collision behaves 
as "repulsive" and the density maximum is likely to be 
located off the trap center resulting in a reduced density. 
In other words, it takes longer time to accumulate high 
enough density for larger <fi. 



A stationary dipolar condensate is known to be stable 
when the atom number is below a critical value. How- 
ever, we find that, even though the total atom number is 
just a fraction of the critical value, during the collision of 
two condensates, the local density fluctuations can still 
induce the collapse. To demonstrate this, we have per- 
formed full three-dimensional numerical simulations for 
typical experimental parameters. We present density dis- 
tributions as function of time for different relative phases 
and compare them with those of two non-interacting con- 
densates. We find that the dipolar interaction modifies 
the density profiles significantly. In addition, we show 
the collapse time as function of the relative phase be- 
tween the two condensates. It turns out that a larger rel- 
ative phase tends to increase the collapse time. We hope 
our study can be helpful to the ongoing experiments with 
degenerate dipolar gases. 




FIG. 6: Linear density distribution p(x) as function of time 
t for <j> = 7r. Blue (red) stands for the minimum (maximum) 
density. 



IV. CONCLUSION 

In conclusion, we have studied the collision of two 
Bose-Einstein condensates with pure dipolar interaction. 
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FIG. 7: Collapse time t c as function of the relative phase <j>. 
Data points with blue square are numerical results. Dashed 
line is a guide to eyes. 



